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Abstract 

A Lagrangean and a set of Feynman rules are presented for Non-Relativistic QCD 
with manifest power counting in the heavy quark velocity v. A regime is identified in 
which energies and momenta are of order Mv. It is neither identical to the ultrasoft 
regime corresponding to radiative processes with energies and momenta of order Mv 2 , 
nor to the potential regime with on shell quarks and gluons providing the Coulomb 
binding. In this soft regime, gluons are on shell, and the quark propagator becomes 
static. Examples show that it contributes to one- and two- loop corrections of scattering 
and production amplitudes near threshold. The results are readily generalised to any 
effective field theory with more than one low energy scale. 
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1 Introduction 



Why a talk on Non-Relativistic QCD [0 || in a workshop dedicated to nuclear physics with 
effective field theories? First, a few words on the background: NRQCD (which in the follow- 
ing also includes NRQED) uses the fact that systems consisting of two or more heavy quarks 
with mass M become very similar to the hydrogen atom: The Coulomb interaction rules the 
level spacing in Charmonium and Bottomium because a s is small enough for perturbative 
calculations. The relative velocity of the quarks in such systems is v ~ a s (Mv) by virtue 
of the virial theorem, where the scale at which the running coupling is to be taken is the 
inverse Bohr radius Mv of the system. Although a s increases with decreasing Q 2 ~ (Mv) 2 , 
a window between the relativistic perturbative and the confinement regime remains in which 
both a s and v is small (for Bottomium, a s (Mv) « 0.2). In the resulting non-relativistic 
and perturbative framework, potentials and wave functions may be used, so that bound 
state physics is easily accounted for, and calculations of production cross sections, hyperfine 
splittings, lifetimes etc. are much facilitated. 

Like in effective nuclear theories, the NRQCD Lagrangean in terms of the heavy quark 
(anti-quark) bi-spinors Q (Q) and gluons (D^ = <9 M + igA^) 

£nrqcd = Q f ( ia ° - 9Ao)Q + ^ Q ] D 2 Q + JjjL QWQ + ... 

+ ^Q^-BQ +... - (CQt*Q) . (Q^CQ) +... (1) 



2M ^ 4M 2 

- j F ^ a + t( Q )f; u d 2 f^ +...+ £ 
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consists of infinitely many terms constrained only by the symmetries of the theory (here: 
gauge invariance) and is non-renormalisable. Predictive power is nonetheless established 
when only a finite number of terms contribute to a given order in the two expansion param- 
eters^, g and v. Excitations with four-momenta bigger than M are integrated out, giving 
rise e.g. to four-point interactions between quarks (d\ ^ 0). So, the ultraviolet physics is en- 
coded in the coefficients q, dj, e^. An advantage of NRQCD is that these can be determined 
by matching NRQCD matrix elements to their QCD counterparts, as both are perturbative 
in the coupling constant. At tree level, the Foldy-Wouthuysen transformation gives Cj = 1, 
and loop corrections are down by powers of g, the most famous example being the coefficient 
for the Pauli term related to the anomalous magnetic moment of the electron, C3 = . .. 

At one loop level, further coefficients enter, and di — — 1. In contradistinction, the free 
parameters of effective nuclear theories have to be determined from data because it is still 
unclear how to match nuclear theory to QCD. 

Another similarity between NRQCD and effective nuclear theory is the existence of three 
scales: Besides the heavy quark mass M, the typical energy and momentum scales in NRQCD 
are the bound state energy Mv 2 and the relative momentum of the two quarks Mv (i.e. the 
inverse size of the bound state) |[[ |2|. In jVjV scattering, the pion production threshold 
a/Matttv and the pion mass supplement the nucleon mass M^. As noted by Luke and 
Manohar 0, this allows one to formally identify a small parameter in nuclear theories, 
~ v. In both theories, the effective Lagrangean does not exhibit the non-relativistic 

1 For clarity, the two will be distinguished in the following, although v ~ a s . 
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expansion parameter explicitly, so that a power counting scheme has to be established which 
determines uniquely which terms in the Lagrangean must be taken into account to render 
consistent calculations and predictive power to a given order in v. It is at this point that 
NRQCD can serve as a "toy model" for effective nuclear theory (although this grossly under- 
states its value): It will establish what the relevant kinematic regimes and infrared variables 
are in a theory with three (or more) separate scales, and it will demonstrate how to count 
powers of v. 

Velocity power counting in NRQCD and identification of the relevant energy and momentum 
regimes has proven more difficult than previously believed. The first attempt by Lepage and 
co-workers 0, fell shot as working only in the Coulomb gauge and as being complicated and 
incomplete. In a recent article, Beneke and Smirnov M pointed out that the much simpler 
velocity rescaling rules proposed by Luke and Manohar for Coulomb interactions ||, and by 
Grinstein and Rothstein for bremsstrahlung processes ||, as united by Luke and Savage ||, 
and by Labelle's power counting scheme in time ordered perturbation theory j7|, do not 
reproduce the correct behaviour of the two gluon exchange contribution to Coulomb scat- 
tering between non-relativistic particles near threshold. This has cast some doubt whether 
NRQCD, especially in its dimensionally regularised version ||, can be formulated using a 
self-consistent low energy Lagrangean. A recent article || has resolved this conflict in a toy 
model, and this contribution presents the extension to NRQCD. 

It is organised as follows: In Sect. |2], the relevant regimes of NRQCD are identified, 
extending the formalism of Luke and Savage || by the soft regime of Beneke and Smirnov J|. 
Sect, ^proposes the rescaling rules necessary for a Lagrangean with manifest velocity power 
counting and gives the vertex and loop velocity power counting rules. An example in Sect. |] 
establishes the necessity of the new, soft regime introduced in Sect. ^ It is essential for the 
correct reproduction of the infrared behaviour of QCD in NRQCD scattering amplitudes. 
Summary and outlook conclude the article. There is also some overlap between the topic of 
this and Mike Luke's talk at this workshop, although I tried to set slightly different priorities. 



2 Idea of Dimensionally Regularised NRQCD 

The NRQCD propagators are 

Q = , A" , (2) 

where T = po — M = ^ + ... is the kinetic energy of the quark. "Num" are numerators 
containing the appropriate colour, Dirac and flavour indices and the gauge fixing term for 
the gluons, all of which are unimportant for the considerations in this section. 

Cuts and poles in scattering amplitudes close to threshold stem from bound states and on- 
shell propagation of particles in intermediate states. They give rise to infrared divergences, 
and in general dominate contributions to scattering amplitudes. With the two scales at 
hand, and energies and momenta being of either scale, three regimes are identified in which 
either the quark or the gluon in (fj) is on shell: 

soft regime: A£ : & ~ \k\ ~ Mv , 
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potential regime: Q p : T ~ Mu 2 , |p| ~ , (3) 
ultrasoft regime: A\^\ k ~ ~ Mv 2 

Ultrasoft gluons -A^ are emitted as bremsstrahlung or from excited states in the bound 
system, and hence physical. Soft gluons A£ do not describe bremsstrahlung: Because in- 
and outgoing quarks Q p are close to their mass shell, they have an energy of order Mv 2 . 
Therefore, overall energy conservation forbids all processes with outgoing soft gluons but 
without ingoing ones, and vice versa, as their energy is of order Mv. 

The list of particles is not yet complete: In a bound system, one needs gluons which 
change the quark momenta but keep them close to their mass shell, relating the (instanta- 
neous) Coulomb interaction: 

A£ : k ~ Mv 2 , \k\ ~ Mv (4) 

So far, only potential gluons and quarks, and ultrasoft gluons had been identified in the 
literature of power counting in NRQCD [R |5|, [JJ. That the soft regime was overlooked 
cast doubts on the completeness of NRQCD after Beneke and Smirnov ||] demonstrated 
its relevance near threshold in explicit one- and two-loop calculations. Here, the fields 
representing a non-relativistic quark or gluon came naturally by identifying all possible 
particle poles in the non-relativistic propagators, given the two scales at hand. 

When a soft gluon A£ couples to a potential quark Q p , the outgoing quark is far off its 
mass shell and carries energy and momentum of order Mv. Therefore, consistency requires 
the existence of quarks in the soft regime as well, 

Q s : T ~ \p\ ~ Mv . (5) 

As the potential quark has a much smaller energy than either of the soft particles, it can - 
by the uncertainty relation - not resolve the precise time at which the soft quark emits or 
absorbs the soft gluon. So, we expect a "temporal" multipole expansion to be associated 
with this vertex. In general, the coupling between particles of different regimes will not be 
point-like but contain multipole expansions for the particle belonging to the weaker kinematic 
regime. For the coupling of potential quarks to ultrasoft gluons, this has been observed by 
Grinstein and Rothstein ||, and by Labelle j7]. 

Propagators will also be different from regime to regime. In order to clarify the relation 
of the presentation here and the work by Beneke and Smirnov [Q, let us consider a typical 
loop integral in NRQCD. One may expand the integrand about the various saddle points, 
i.e. about the values of the loop-momentum q where particles become on shell. For example, 
expanding about a saddle point coming from a physical gluon at the soft scale, a quark 
propagator may be expanded as (T p ~ ~ Mv 2 <C qo ~ \q\ ~ Mv) 

i i_ i • frp (P + Q) 2 \ i rfi x 

So, Q s is expected to become static to lowest order, and the higher order terms in the 
expansion can be interpreted as insertions into the soft quark propagator, or (as mentioned 
above and to be confirmed below) as resulting from an energy multipole expansion which 
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modifies the vertex rules. As the energy of potential gluons is much smaller than their 
momentum, the A p -propagator is expected to become instantaneous for similar reasons. 

With these five fields Q s , Q p , A£ , A^ , A^ representing quarks and gluons in the three dif- 
ferent non-relativistic regimes, soft, potential and ultrasoft, NRQCD becomes self-consistent. 
An ultrasoft quark (which would have a static propagator) is not relevant for this paper. It 
is hence not considered, as is a fourth ( "exceptional" ) regime in which momenta are of the 
order Mv 2 and energies of the order Mv or any regime in which one of the scales is set by 
M. They do not derive from poles in propagators, and hence will be relevant only under 
"exceptional" circumstances. A future publication || has to prove that the particle content 
outlined is not only consistent but complete. 

It is worth noticing that the particles of the soft regime can neither be mimicked by 
potential gluon exchange, nor by contact terms generated by integrating out the ultraviolet 
modes: Fields in the soft regime have momenta of the same order as the momenta of the 
potential regime, but much higher energies. Therefore, seen from the potential scale they 
describe instantaneous but non-local interactions, as pointed out by Beneke and Smirnov J|. 
Integrating out the scale Mv, one arrives at soft gluons and quarks as point-like multi-quark 
interactions in the ultrasoft regime. The physics of potential quarks and gluons will still 
have to be described by spatially local, but non- instantaneous interactions. This suggests 
once more that there is no overlap between interactions and particles in different regimes. 

Finally, the regularisation scheme must be chosen such that the three kinematic regimes 
still do not overlap, i.e. such that expansion around one saddle point in the loop integral 
does not obtain any contribution from other saddle points and regimes. One might use an 
energy and momentum cutoff separating the soft from the potential, and the potential from 
the ultrasoft regime, but the integrals encountered can in general not be performed analyti- 
cally. Furthermore, cutoff regularisation usually jeopardises power counting and symmetries, 
and introduces unphysical power divergences as the (unphysical) cutoff is removed. In con- 
tradistinction, using dimensional regularisation after the saddle point expansion preserves 
power counting and gauge symmetry. Its homogeneity @ guarantees that contributions from 
different saddle points and regimes do not overlap (A simple example can be found in Ref. 
||.). Therefore, dimensional regularisation will be the method of choice in the example of 
Sect. |. 



3 Velocity Power Counting 
3.1 Rescaling Rules and Propagators 

In order to establish explicit velocity power counting in the NRQCD Lagrangean, one rescales 
the space-time coordinates such that typical momenta in either regime are dimensionless, 
as proposed by Luke and Manohar || for the potential regime, and by Grinstein and Roth- 
stein [H for the ultrasoft one: 

soft: t = (Mvy 1 T s , x = (Mvy 1 X s , 
potential: t = (Mv 2 )' 1 T u , x = (Mv)' 1 X s , (7) 
ultrasoft: t = (Mv 2 )' 1 T u , x = {Mv 2 )- 1 X u . 
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For the propagator terms in the NRQCD Lagrangean to be properly normalised, one sets 
for the representatives of the gluons in the three regimes 



soft 
potential 
ultrasoft 

and for the quark representatives 

soft: 
potential: 



A2(x,t) = (Mv) A»(X S ,T S ) 
A»(x,t) = (Mvi) A»(X 8 ,T U ) 
A£(x,t) = (Mv 2 ) A£(X U ,T U ) 

Q s (x,t) = (Mv)l Q S (X S ,T S ) 
Q p (x,t) = (Mv)% Q P (X S ,T U ) 



(8) 



(9) 



The rescaled free quark Lagrangean reads then 

soft: d 3 x s dr s g s t(id + ^)g s , 

1 



potential: d 3 X s dT u Q p f (id + - d 2 ) Q 



(10) 
fill 



Here, as in the following, the positions of the fields have been left out whenever they coincide 
with the variables of the volume element. Derivatives are to be taken with respect to the 
rescaled variables of the volume element. 

The gauge fixing term was included in the NRQCD Lagrangean (|IJ) because the decom- 
position of the Lagrangean into a free and an interaction part is gauge dependent. Usually, 
the Coulomb gauge d ■ A = is chosen in NRQCD, but Luke and Savage |3j showed how to 
establish explicit velocity power counting in any gauge. Because of the difference between 
canonical and physical momentum, it is important to specify the gauge before identifying to 
which order in v a certain regime in the Lagrangean contributes, as seen shortly. Still, the 
classification of the three kinematic regimes (|3|) itself relied only on the typical excitation 
energy and momentum, and hence on gauge invariant quantities, and the denominator in 
the gluon propagator (H) is gauge independent. 

The rescaled free gluon Lagrangean in the Lorentz gauge reads for example 



soft: 
potential: 

ultrasoft: 



d s X s dT s - A% 



a 



a: 



d 3 X s dT u -A^[g^(v 2 d 2 -d 2 ) 



d 3 X u dT u -A^[d 2 gflI/ 



){v8^ 9o + 5^di)(v5 u0 do + 5 ui di)\A^ 

l--)d^d v }A u u , 

a J 



(colour indices suppressed), while in the Coulomb gauge 



soft: 
potential: 
ultrasoft: 



2 

d 3 X s dT u i 



d 6 X s dT s ^ Ai, s - d 2 )5 tJ - didj A j)t 



A, P d 2 A, P + Av( d 5 H - didj - v'd^Aj 



d 3 X s dT s - Ai >u [{^ - d>)5ii - (h<h Aj,; 



(12) 
(13) 

(14) 

(15) 
(16) 
(17) 
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The (un-rescaled) Coulomb gauge propagators are therefore (Dirac and colour indices sup- 
pressed, 5% = <ytf - 

(t,p) i t k j isit 



tr 



T + ie ' h ' A; 2 + ie 

i 



(19) 



—A; 2 + ie 1 —A; 2 + ie 

1 °tr 



k 2 + ie 



(20) 



As expected, the soft quark becomes static and the potential gluon becomes instantaneous in 
both gauges. In order to maintain velocity power counting, corrections of order v or higher 
must be treated as insertions, represented in Coulomb gauge by the (un-rescaled) Feynman 
rules 

(T,p) _ 2 k,A 

)( = = 0{v l ) , 1— X— ' = +ik 2 5 l3 = 0(v 2 ) . (21) 

The Lorentz gauge propagators and insertions are written down straightforwardly, too ||, 
and look very similar to the Coulomb gauge result, especially for a — 1. As seen from 
(|l5Hi7l) ; the choice of the Coulomb gauge makes A instantaneous, and hence it contributes 



in the potential regime, only. Since in this gauge, A solely mediates the instantaneous 
Coulomb potential (physical fields are transverse by virtue of GauB' law), this result was to 
be expected. The field A p is associated with retardation effects like spin-orbit coupling and 
the Darwin term in (|I|). The advantages of having couplings between Aq and the other fields 
only for potential Aq and of having no insertions in the Aq propagator is balanced by the 
demand for non-multiplicative renormalisation of the Coulomb gauge. The Lorentz gauge 
may hence facilitate some calculations although the number of diagrams is larger, as will be 
seen shortly. 

Except for the physical gluons A£ and A^ there is no distinction between Feynman and 
retarded propagators in NRQCD: Antiparticle propagation has been eliminated by the field 
transformation from the relativistic to the non-relativistic Lagrangean, and both propagators 
have maximal support for on-shell particles, the Feynman propagator outside the light cone 
vanishing like e~ M . Feynman's perturbation theory becomes more convenient than the time- 
ordered formalism, as less diagrams have to be calculated. 



3.2 Vertex Rules 

By experience, particles in the various regimes couple: On-shell (potential) quarks radiate 
bremsstrahlung (ultrasoft) gluons. In general, one must allow all couplings between the 
various regimes which obey "scale conservation" for both energies and momenta. They must 
be conserved within each regime to the order in v one works. This will exclude for example 
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the coupling of two potential quarks (T ~ Mv 2 ) to one soft gluon (q ~ Mv), but not to 
two soft gluons via the Q^A ■ AQ term of the Lagrangean ([!]). 

As an example, consider a bremsstrahlung-like process: the radiation of a soft gluon off 
a soft quark, resulting in a potential quark. The rescaled interaction Lagrangean reads for 
the vector coupling 



d 3 X dT s 



-igvQ s \X s ,T s )d-A s (X s ,T s )Q p (X s ,vT s ) . (22) 



Note that the scaling regime of the volume element is set by the particle with the highest 
momentum and energy. The Feynman rule for this vertex is hence 

(T,p) (T>,p>) 

-ig(p-p')i(2n) 4 5^(p+p' + q)x 



q,i 



exp (T' 8(T + q )} = 0(v e 



One sees that technically, the energy multipole expansion expected in Sect. |2] comes from the 
different scaling of x and t in the three regimes. The factor e^ symbolises that the multipole 
expansion corresponds term by term to an expansion in v. It should be truncated at the 
desired order in v. 

Amongst the fields introduced, six interactions are allowed within and between the var- 
ious regimes for the vector coupling, and two (six) for the scalar coupling in the Coulomb 
(Lorentz) gauge. Their v counting is presented in tables [I] and |^. Note that - albeit both 
describing interactions with physical gluons - soft and ultrasoft couplings occur at different 
orders in v, and obey different multipole expansion rules. On the level of the vertex rules, 
double counting is prevented by the fact that in addition to most of the propagators, all 
vertices are distinct because of different multipole expansions. Velocity power counting for 
other vertices is obtained again by rescaling and multipole expansion. For example, the rules 
for the Fermi term £ int = ^ Q^a ■ BQ are identical to those of table [TJ. 



Table 1: Velocity power counting and vertices for the interaction Lagrangean — 1| Q^d ■ AQ 
in both Lorentz and Coulomb gauge. 

















Vertex 


I 

1 A 
1 






% 


1 A 
1 




v power 




V 


V 


V 


3 

V 2 


v 2 



Table 2: Velocity power counting and vertices for the interaction Lagrangean —gQ^A Q in 
the Lorentz gauge. In the Coulomb gauge, only the hrst two diagrams exist. 

















Vertex 


I 

| A 
1 


| A 
1 










v power 


1 




V° 


V° 


v° 


V 



Using the equations of motion, a temporal multipole expansion may be re-written such 
that the energy becomes conserved at the vertex. Now, both soft and potential or ultrasoft 



s 



energies are present in the propagators, making it necessary to expand it in ultrasoft and 
potential energies. An example would be to restate the vertex (^) as 

| = -ig(p-p')d^yS(T + T + q )5^(p+p' + q) , (24) 

and the soft propagator to contain insertions 0(v) for potential energies T" 

(-T=qp+T',p) j /_T' 

qo + ie ^ V Qo 

This gives us the first part of the correspondence to Beneke and Smirnov's threshold expan- 
sion in the example propagator (^|), the second one coming from the soft quark insertion 
(|2"1"D. The same can be shown for the momentum-non-conserving vertices, too. 

It is also interesting to note that there is no choice but to assign one and the same coupling 
strength g to each interaction. Different couplings for one vertex in different regimes are not 
allowed. This is to be expected, as the connection to Beneke and Smirnov's threshold 
expansion |3| demonstrated that the fields in the various regimes are representatives of one 
and the same non-relativistic particle, whose interactions are fixed by the non-relativistic 
Lagrangean (|26|) . 

In the renormalisation group approach, there is only one relevant coupling (i.e. only one 
which dominates at zero velocity): As expected, it is the Q p Q p A p0 coupling providing the 
binding. In the Coulomb gauge, all other couplings and insertions are irrelevant, while the 
Lorentz gauge exhibits three marginal couplings: Q P Q P A U ^ , Q s Q s A S fi and Q s Q P A S fi. 




3.3 Loop Rules 

The velocity power counting is not yet complete. As one sees from the volume element 
used in (0), the vertex rules for the soft regime count powers of v with respect to the soft 
regime. One hence retrieves the velocity power counting of Heavy Quark Effective Theory ]10, 



|TT|| (HQET), in which the interactions between one heavy (and hence static) and one or 
several light quarks are describedP]. HQET becomes a sub-set of NRQCD, complemented by 
interactions between soft (HQET) and potential or ultrasoft particles. 

In NRQCD with two potential quarks as initial and final states, the soft regime can occur 
only inside loops, as noted above. Therefore, the power counting in the soft sub-graph has 
to be transfered to the potential regime. Because soft loop momenta scale like [d 4 /c s ] ~ v 4 , 
while potential ones like [d 4 A;p] ~ v 5 , each largest sub-graph which contains only soft quarks 
and no potential ones (a "soft blob") is enhanced by an additional factor -. 

As an example, consider the graphs of Fig. [I]: Using the Lorentz gauge, vertex power 
counting gives that the leading contribution is from the exchange of two potential gluons, 
coupled via Q^QA Q . There are four such vertices, so the diagram is 0(g 4 v~ 2 ) (table |2|). The 
next two diagrams are O(g 4 v ) and O(g 6 v ) from the vertex power counting, but another 
factor - must be included because there is one soft blob in the diagrams. The intermediate 
couplings in the third diagram take place in the soft regime and hence are counted in that 



2 Usually, HQET counts inverse powers of mass in the Lagrangean, but because in the soft regime Mv 
const., the two approaches are actually equivalent. 
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regime. The last diagram, in which two soft blobs are separated by the propagation of two 
potential quarks is O(g 8 v ) from the vertices, and the loop counting gives a factor Each 
soft blob contributes at least four orders of g, but only one inverse power of v ~ g 2 . Power 
counting is preserved. These velocity power counting rules in loops are verified in explicit 
calculations of the exemplary graphs (see also below), but a rigorous derivation is left for a 
future publication [|]. 




Figure 1: Power counting with soft loops. The loops in the second and third diagram obtain 
an inverse power of v , the last diagram of v 2 in addition to the power counting following 
from the vertex rules (fj). 

With rescaling, multipole expansion and loop counting, the velocity power counting rules 
are established, and one can now proceed to check the validity of the proposed Lagrangean by 
matching NRQCD to the relativistic theory in an example given by Beneke and Smirnov M. 

4 A Model Calculation 

By construction, NRQCD and QCD must agree in the infrared limit, and especially in the 
structure of collinear (infrared) divergences. Matching NRQCD to the low velocity limit of 
QCD will therefore confirm that the power counting proposed is correct and that the soft 
regime is relevant. In QCD, the lowest order infrared divergent contribution in the Coulomb 




Figure 2: Leading soft (left) and potential (right) planar contribution to quark-quark scat- 
tering in the Coulomb gauge. 

gauge which comes from the soft regime is the one loop graph Fig. |2| of order g 4 v, while the 
leading contribution at this order in g comes from the two Coulomb gluon exchange and is 
order v~ 2 . The leading diagram mediating retardation effects (A p couplings) in this order in 
g is 0(v 2 ). At high enough order in g, soft physical gluon exchange will therefore be more 
significant than retardation effects. The number of diagrams to be considered in the case at 
hand is large enough^ to make one long for a simpler example. 

3 Actually, most of them are zero, but to establish a pattern which convinces one of that is not the topic 
of this contribution || . 
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For the sake of simplicity, let us - following Beneke and Smirnov || - deal with a toy 
model NRQFT Lagrangean 



£nrqft = &(id + — -g Cl A)<S> + -{d,A){d»A) + c 2 (V$) + . . . ( 26 ) 

of a heavy, complex scalar field $ with mass M coupled to a massless, real scalar A. The 
coupling constant g has been chosen dimensionless. In a slight abuse of language, $ will still 
be referred to as "quark" and A as "gluon". The coefficients are again to be determined 
by matching relativistic and non-relativistic scattering amplitudes. This Lagrangean is very 
similar to the A -part of the NRQCD Lagrangean in Lorentz gauge. Especially, the vertex 
and loop power counting is identical to the one of table ||. 

The first soft non-zero contribution in this toy model comes from the two gluon direct 
exchange diagram of Fig. |3] calculated by Beneke and Smirnov Q using threshold expansion. 
The Mandelstam variable t = — (p — p') 2 describes the momentum transfer in the center of 
mass system, y = — (p) 2 oc — v 2 the relative four-momentum squared of the ingoing quarks 
as indicator for the thresholdness of the process. The ultraviolet behaviour of this graph 



(TjP) 




(-fc ,-p-fe) 

Figure 3: Planar 0(g 4 ) contributions to Coulomb scattering in the toy model. The four-point 
interaction and insertion diagrams are not displayed. 



is mimicked in NRQFT by the four-fermion exchange with ic 2 = 24 ~2 9 M2 = O(t ,y°), which 
using the rescaling rules is seen to be Oiy 1 ). 

The v4 u v4 u -diagram is of order e v (see table |2|) with a leading loop integral contribution 
(similar to Beneke and Smirnov's Q fl. (32)) 

r _^ i L_ i i (27) 

The diagram is expected to be zero since the ultrasoft gluons do not change the quark 
momenta and therefore the scattering takes place only in the forward direction, p = p'. As 
no scale is present, it indeed vanishes upon employing the on-shell condition for potential 
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quarks, T = ^ to leading orderfj. The A u A p and A p A u contributions (0(^e v )) are zero for 
the same reason. The lowest order contribution to the A p t4 p graph (O(t)) is: 

d* i i i i 



(2tt)" fc2 -ie (p fc) 2 -KT + ko- ffifil + ie T — fco — M + if 



In the light of the discussion at the end of Sect. it is most consistent to perform the 
k integration by dimensional regularisation, using |L2| / j0p = J J^p- ^t^-* > a ~ > 



Split dimensional regularisation was introduced by Leibbrandt and Williams |L3| to cure the 
problems arising from pinch singularities in non-covariant gauges. Here, it has the same 
effect as closing the fco-contour and picking up the quark propagator poles prior to using 
dimensional regularisation in d — 1 Euclidean dimensions. Considering also one insertion 
( |2l"D at the potential gluon lines to achieve 0(v *) accuracy, the result, 

i M + T f 2 , -t n 

7E -ln— ) , (29) 



8-jrt y/y H - d A-kij? 

agrees with Beneke and Smirnov's fl. (31) when one keeps in mind that non-relativistic 
external lines were normalised differently, and that different conventions for dimensionally 
regularised integrals were chosen. Near threshold, the scale is set by the total threshold 
energy 4yr/i 2 = 4(M + T) 2 . 

The soft gluon part is to lowest order {0{v~ 1 )) given by 



/ 



■ n 1 iii 



(2n) d k 2 -k 2 + ie k% - (p - f + W + ie + ie ~k + ie 



which corresponds to Beneke and Smirnov's M fl. (33). Now, split dimensional regularisation 
must be used if no ad-hoc prescription for the pinch singularity at fco = is to be invoked. 
That the pinch is accounted for by potential gluon exchange and hence must be discarded, 
agrees with the intuitive argument that zero four-momentum scattering in QED is mediated 
by a potential only, and no retardation or radiation effects occur. On the other hand, the 
model Lagrangean contains three marginal couplings as seen at the end of Sect. |3]2|, which 
may give finite contributions as energies and momenta of the scattered particles go to zero. 
The result to 0(v l ) exhibits another collinear divergence, 



\ / 'l — t \ l r 'ly ( m z — t \i 

^(l37/-^- ln w) + 5I^i 1 + T(l^-^- ln w)J ■ ( 31 ) 



2y r 2 -t 



AirH M — d ' 4yr/i 2 / 24tt 2 M 2 L ' t ^A-d ,Sh 4?r^ 

and agrees with fl. (36) given by Beneke and Smirnov ||]. The second term comes from 
insertions and multipole expansions to achieve 0(v l ) accuracy. 

It is easy to see that the power counting proposed works. As expected, the potential 
diagram is ^/y oc v stronger that the leading soft contribution, and tyfy oc v 3 stronger than 
the four-fermion interaction. 

In conclusion, the proposed NRQCD power counting and Lagrangean with three different 
kinematic regimes @ reproduces the collinear divergences of the planar graph of the rela- 
tivistic theory only if the soft gluon and the soft quark are accounted for: The four-fermion 

4 Since T — ~ Mv 4 <C k ~ Mv (Mv 2 ) in the potential regime, this is a legitimate expansion. 
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contact interaction produces just a -pa-term, graphs containing ultrasoft gluons were ab- 
sent, and the potential gluon ( P5| ) gave no O(y ) contribution. The coupling strength of the 
$ s A;$p vertex is also seen to be identical to the other vertex coupling strengths, g. 



5 Conclusions and Outlook 



The objective of this contribution was a presentation of the ideas behind explicit velocity 
power counting in dimensionally regularised NRQCD. The identification of three different 
regimes of scale for on-shell particles in NRQCD leads in a natural way to the existence 
of a new quark field and a new gluon field in the soft scaling regime E ~ \p\ ~ Mv. 
In it, quarks are static and gluons on shell, and HQET becomes a sub-set of NRQCD. 
Neither of the five fields in the three regimes should be thought of as "physical particles". 
Rather, they represent the "true" quark and gluon in the respective regimes as the infrared- 
relevant degrees of freedom. None of the regimes overlap. An NRQCD Lagrangean has been 
proposed which leads to the correct behaviour of scattering and production amplitudes. It 
establishes explicit velocity power counting which is preserved to all orders in perturbation 
theory, once dimensional regularisation is chosen to complete the theory. The reason is non- 
commutativity of the expansion in small parameters with dimensionally regularised integrals. 

I would like to stress that the diagrammatic threshold expansion derived here allows for 
a more automatic and intuitive approach and makes it easier to determine the order in v to 
which a certain graph contributes than Beneke and Smirnov's way Also, the NRQCD 
Lagrangean can easily be applied to bound state problems. An investigation of the influence 
of soft quarks and gluons on bound state calculations in NRQED and NRQCD is important 
because - as seen at the beginning of Sect. ^, their contribution at 0{g i ) and higher becomes 
stronger than retardation effects. As the threshold expansion of Beneke and Smirnov starts 
in a relativistic setting, it may formally be harder to treat bound states there. 

Coming back to the topic of this workshop, effective nuclear theories, NRQCD shows 
how to establish a power counting in any effective field theory with several low energy scales: 
First, identify the combinations of scales in which particles become on shell by looking at 
the denominators of the various propagators. This gives the scaling regimes. Then, the 
Lagrangean is rescaled to dimensionless fields in each regime to exhibit the vertex and loop 
power counting rules. A priori, all couplings obeying scale conservation are allowed. 

The problem with effective NN scattering is not that the three scales M N , y/Mjqm^ , 



are separated only by powers of yj%£ « 0.4. Indeed, Kaplan, Savage and Wise |14L na 
obtain very promising results for scales much smaller than the pion mass (See also David 
Kaplan's talk in this workshop, and Martin Savage's contribution on the inclusion of what 
in the language of this article would be ultrasoft pions.). The difficulty is that the scale 
^/M^m^ pa 360MeV at which the soft regime becomes relevant is larger than the NN 
scattering expansion parameter Ann ~ 300MeV. How to overcome this is another interesting 
topic for the future. 
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